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PREFACE 


Tlie  investigation  reported  herein  was  conducted  hy  personnel  of 
the  Soil  ]>ynaxttics  Division  (SDD',  Soils  and  Pavements  Laboratory  (S&PL), 
U.  S.  Army  Engineer  Wate-ways  Experiment  Station  (WES),  as  a Tjart  of  DA 
Project  No.  UA16i10?B52E,  "Fragment  and  Px'ojectile  Penetration  Resistance 
of  Soils,"  Task  0^4,  Work  Unit  13- 

The  theoretical  developments  weie  originated  by  Mr.  R.  S.  Bernard 
during  the  period  July  19T^  to  January  19T5  under  the  supervision  of 
Dr.  J.  G.  Jackson,  Jr.,  Chief,  SDD.  Messrs.  J.  P.  Sale  and  R.  G.  Ahlvin 
were  Chief  and  Assistant  Chief,  S&PL,  respectively.  SP5  D.  C.  Creighton 
developed  and  executed  the  necessary  computer  codes,  ana  Drs.  B.  Rohani 
and  P.  Fo  Hadala  provided  continuous  technical  and  moral  support.  This 
report  was  prepared  by  Mr.  Bernard. 

COL  G.  H.  Kilt,  CE,  was  Director  cf  WES  during  the  investigation 
and  at  the  time  of  publication  of  this  report.  Mr.  F.  R.  Brown  was 
Technical  Director. 
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Note: 


1 kilobar  = IjOOO  bars  = lUj500  poimds  per  square  inch  = 100 
megapascals » 
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DEVELOmfflT  OF  A PHOJECTILE  PMETMTIOl  THEORY 


DES*  PMETRATIOI  THEORY  FOR  HOMOGEHEOUS  AID  LAYERED  TARGETS 

PART  I ■.  INTRODUCTION 
Background 


1.  Projectile  penetration  has  been  the  object  of  theoretical  and 
e»5)eriiaental  investigations  for  at  least  two  centtiries.  However ^ until 
recent  years,  reliable  penetration  predictions  were  achieved  only  by 
means  of  empirical  equations  using  various  target  resistance  functions 
and  pseudoconstants  which  were  drawn  diiectly  from  penetration  experi- 
ment-^, The  first  successful  analysis  on  purely  theoretical  groiinds  was 
appax^ntly  due  to  Bishop,  Hill,  and  Mott,^  who  related  quasi-static 
punch  indentation  to  the  static  expansion  of  spherical  and  cylindrical 
cavities  in  an  infinite  medium.  In  a later  development  which  incorp^ora- 

ted  an  ^ hoc  analogy  with  the  dynamic  expansion  of  a spherical  cavity 

2 3 

in  an  infinite  medium,  Goodier  formulated  a dynamic  penetration  theory 

for  rigid  spherical  projectiles  which  is  applicable  to  elastic-plastic 

strain-hardening  targets.  In  subsequent  efforts,  attempts  were  made  to 

U 

account  for  the  effects  of  target  compressibility  and  projectile  nose 
shape, ^ but  the  substance  of  the  Goodier  penetration  theory  remained 
unchanged. 

2,  A useful  penetration  theory  must  be  fairly  simple  in  order 
that  reliable  predictions  can  be  made  with  a minimal  expenditure  of  time 
and  money.  A number  of  simplifying  assumptions  must  be  made  in  the 
development  of  such  a theory^  and  since  the  objective  is  simply  to  ob- 
tain an  equation  of  motion  for  the  projectile,  only  the  gross  features 

of  the  target  response  may  be  taken  into  account.  An  extreme  alternative 
to  the  simple  theoretical  approach  is  provided  by  two-dimensional  (2D) 
finite-difference  analyses  which  generally  require  large  expend! t\ires  in 
man-hours  and  computer  time.  IT.-v-'rtheless , a 2D  finite-difference 
solution  has  the  advantage  that  the  details  of  projectile  and  target 


It 


betiavior  can  be  closely  examined;  and  while  an  approach  of  this  kind  niay 
be  too  costly  and  cvrabersome  for  large-scale  parameter  studies,  it  can 
be  extremely  valuable  for  assessing  the  relative  effects  of  various 
mechanisms  in  the  penetiation  process. 

3.  At  the  present  time,  there  is  widespread  need  for  a ti’actable 
theory  which  can  be  used  to  study  the  penetration  of  soil,  rock,  con- 
crete, and  metal  targets  by  projectiles  such  as  bombs,  bullets,  shell 
fra@nents,  and  missiles.  Such  a theory  is  needed  from  both  the  offen- 
sive and  defensive  viewpoints  since  a theory  which  analyzes  the  pene- 
tration process  can  also  be  used  to  determine  the  conditions  under  which 
final  penetration  depth  is  maximized  or  minimized,  respectively.  In 
July  1973,  a small  research  effort  was  initiated  under  the  sponsorship 
of  the  Office,  Chief  of  Engineers  (OCE),  to  develop  a penetration  theory 
which  would  account  for  the  penetration  of  homogeneous  and  layered  tar- 
gets by  rigid  axisymmetric  projectiles.  In  the  period  Jiily  1973  to  June 
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197^,  the  previously  existing  Goodier  theory  for  spherical  projectiles  ’ 
was  extended  for  projectiles  with  conical  and  ogival  nose  shapes  at 
shallow  to  moderate  depths,  as  reported  in  Reference  5*  In  the  present 
report  this  theory  is  modified  and  extended  to  account  for  deep  penetra- 
tion of  homogeneous  and  layered  targets  at  very  high  impact  velocities. 
However,  it  is  emphasized  that  the  modified  theory  is  still  applicable 
to  shallow  penetration  problems  and  produces  approximately  the  same 
results  as  the  unmodified  theory  at  low  impact  velocities.  The  theory 
developed  herein  supersedes  all  preceding  work  and  is  applicable  for 
shallow  and  deep  penetration  in  cohesive  targets.  The  developments 
presented  in  this  report  are  the  results  of  efforts  conducted  luider  the 
OCE  program  dioring  the  period  July  197^  to  January  1975. 


Purpose 


U.  The  pxirpose  of  this  investigation  is  to  build  upon  existing 
penetration  theory^  to  develop  a more  general  theory  for  cohesive  targets 
which  is  valid  for  deep  penetration  at  very  high  impact  velocities  as 
well  as  for  shallow  penetration  at  low  impact  velocities.  As  in 


5 


previous  i#ork,  the  fundamental  baais  for  the  penetration  theory  is  the 

* 3„5 

dynamie  cavity  expansion  theoiy,  and  the  modification  for  high=speed, 
deep  penetration  xs  accompli through  a po&txilated  rels.tion  between 
target  penetration  resistance  and  the  ”sclid  Reynolds  number,"  which  was 
first  introduced  in  Reference  5- 


Scope 

5.  This  report  is  intended  to  be  e.  self-contained  document  on  the 
subject  of  projectile  penetration  theory  for  cohesive  targets,  and  it  is 
structured  such  that  each  part  may  be  read  independently  or  in  sequence, 
depending  on  the  interests  of  the  reader.  The  cavity  expansion  theory 
is  discussed  in  Part  II.  Previous  work  in  shallow  penetration  theory  is 
outlined  in  Part  III,  followed  by  the  modifications  whei'eby  the  theory 
is  extended  for  deep  penetration.  Practical  applications  of  the  theory 
are  presented  in  Part  IV,  and  conclusions  and  recommendations  are  stated 
in  Part  V. 


PART  IJ:  CAVITY  EXPAifSION  THEORY 


Background 


6,  The  dynamic  expansion  of  spherical  cavities  in  solids  has  been 

2-5  3 

investigated  by  several  authors.  Goodier  used  theory  of  dy- 
namic  cavity  expansion  in  an  incompressible  strain~>hardening  material'^ 
as  the  basis  of  a penetration  theory  for  rigid  spherical  projectiles  and 
homogeneous  targets.  In  Reference  5»  the  cavity  expansion  theory  was 
extended  to  concentrically  layered  compressible  media  and  was  used  as 
the  basis  of  a projectile  penetration  theory  for  layered  targets.  The 
cavity  expansion  theory  developed  in  Reference  5 represents  a first 
attempt  at  a theory  for  concentrically  layered  media  and  is  more  compli- 
cated than  necessary.  In  the  following  discussion,  a simplified  version 
of  the  cavity  exansion  theory  for  concentrically  layered  media  will  be 

presented  and  w provide  a smoother  transition  to  the  penetration 

theory  than  does  the  more  complicated  version. 


Problem  Formiolation 


7.  Consider  a spherical  cavity  of  radius  a(t)  surrounded  by 
two  concentric  layers  of  different  materials,  as  shown  in  Figure  1. 
These  materials  exhibit  elastic-plastic  response  with  linear  strain- 
hardening under  shear  stress,  as  illustrated  in  Figure  2.  Both 


Figure  1.  Spherical  cavity 
expansion  in  a concentrically 
layered  medium 
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Figure  2.  Material  behavior  in  shear 

1+ 

materials  exhibit  ideal  locking  compressibility  under  hydrostatic 
stress  such  that  transition  from  the  elastic  to  plastic  states  is  ac- 
companied by  a small  but  finite  volumetric  strain  known  as  the  locking 
strain 


e 


{!) 


where  and  are  the  material  densities  in  the  elastic  and  locked 

plastic  states,  respectively,  and  |e^|  <<  1 . Otherwise,  the  materials 
are  incompressible,  and  the  regions  of  elastic  and  locked  plastic  be- 
havior are  separated  by  a weak  plastic  shock  front  of  radius  b(t)  , 
where  t denotes  time.  Accordingly,  the  locking  condition  on  strain  is 


, a < r < b 
f 0 , r > b 


(2) 


where  e^(r,  t)  and  £0(^»  t)  are  normal  strains  in  the  radial  and  cir- 
cumferential directions,  respectively,  with  r as  the  radial  coordinate. 
The  associated  strain  rates  are  related  to  the  outward  radial  particle 
velocity  v(r,  t)  by 
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® V 
6.  a — 

0 r 


(3) 


e 
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8r 


where  a dot  above  any  quantity  denotes  a time  derivative » Equations  2 
and  3 give  rise  to 


v(r,  t ) 


a < r < "b 


r > b 


ih) 


where  fn(t)  and  f (t)  are  unknown  functions.  The  stress=-strain 
£ e ^ 

relation  for  the  material  in  shear  is 


- 3 ^ ^ 'b 

(5) 

2ESq  5 r > b 

where  cr^(r5  t)  and  t)  are  normal  stresses  in  the  radial  and 

circumferential  directions  and  Y , Eg  and  E_^  are  the  yield  strength, 
elastic  modulus,  and  strain- hardening  modulus,  respectively,  for  the 
material  in  compression.  The  geometry  is  spherically  symmetric,  and  the 
material  equation  of  motion  is 


o 

r ^ 2 f 

- — + — (a 
8r  r r 


■Jg)  = P 


/3v  . 3v 
\3t  ^ 


3v\ 

3r/ 


(6) 


where  p is  the  local  density.  After  combining  Equations  U,  5»  and  6, 
the  resulting  equation  is  integrated  from  r = “>  to  r = a(t)  , subject 
to  the  condition  that  the  material  is  stress-free  as  r ->  <»  . The  equa- 
tion for  the  compressive  normal  stress  p = ‘-cr^,[a(t)]  at  the  cavity 
surface  is  then  found  to  be; 
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+ 2 I <Y  + 2E, 


CO 


where  Aa,  is  the  discontinuity  in  a at  r = b(t),  given  by 

D I* 


_ J^)  J-) 

- % - % 

and  and  ^ are  the  lin»iting  values  of  o as  r 

D b r 

the  elastic  and  locked  plastic  regions,  respectively. 


Problem  Solution 

8.  In  order  to  evaluate  Equation  T and  thereby  obtain  an  explicit 
relation  between  p and  a(t)  , relations  must  be  established  between 
a(t)  and  fg^(t)  , i'g('t)  > , h(t)  , b(t),  and  £0(r,  t)  . 

The  first  of  Equations  4 is  evaluated  at  r = a leading  to 

fp  = a"a  (9) 

The  jump  condition  for  conservation  of  mass  is  applied  at  r = b so  that 

p^[b  - v^""^]  = pQ[b  ~ v^'*’^]  (10) 

( + ) (“) 

where  v and  v are  the  limiting  values  of  v as  r b in  the 
elastic  and  locked  plastic  regions,  respectively.  After  evaluating 
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Equation  ^ at  r - b and  combining  the  re-isults  with  Equations  9 and  10 
the  following  is  obtained: 


2®  ?® 
a a -•  ab  "b 

1 - a 


(11) 


where  a is  the  material  compaction  coefficient  defined  by 


a 


(12) 


which  is  limited  to  a <<  1 since  attention  has  been  restricted  to 
|e^|  <<  1 = The  jump  condition  for  conservation  of  momentum  is  applied 
at  r = b so  that 


PqV 


(^) 


(13) 


Evaluating  Equation  4 at  r = b and  then  combining  the  results  with 
Equations  9®  11®  12 9 and  13  leads  to 


®2/-  2/^2x2 

apg^b  (1  - a /b  ) 

1 - a 


(lU) 


For  large  deformations,  the  circtitnferential  strain  can  be  appropriately 
expressed  as 


= In  7-  (15) 

o 

where  is  the  initial  position  of  a particle  located  at  position  r . 

Restricting  attention  to  situations  wheie  Y <<  E , E^  <<  E , and 
le^j  « 3Y/2E^  , then 


EglbCt)]  = In  s 3 (16) 

c 

subject  to  P <<  1 where  b^  is  the  initial  position  of  a particle 
located  at  r = b , and 
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0 = 2 (e  - 


9.  Attention  is  now  restricted  to  materials  which  have  identical 
values  of  a and  . (This  restriction  can  he  relaxed  hy  establishing 
appropriate  composite  values  of  a and  for  the  two  layers.) 

Approximations  for  t)  » h(t)  , and  h(t)  will  now  he  obtained 

subject  to  the  following  assumptions: 

a.  Prior  to  the  arrival  of  the  plastic  front  at  the  initial 

layer  interface  location  (i.e.,  h < h ),  the  entire  second 
layer  is  in  the  elastic  state,  and  the'^location  of  the 
interface  is  approximately  h(t)  — • 

After  the  plastic  front  arrives  at  the  initial  layer  in- 
terface position  (i.e,,  b ^h^),  the  locked  plastic  region 
extends  into  the  second  layer,  and  the  motion  of  the  in- 
terface is  no  longer  neglible. 

10.  Considering  first  the  situation  when  b < h and 
a^(t)  >>  a^(0)  , conservation  of  mass  in  the  locked  plastic  regior^ 
gives  rise  to 


3 3 3 

r - a a (1  - a)r^  , a < r < b 


(18) 


Rearranging  Equation  l8  and  then  combining  the  result  with  Equation  15 
yields  the  expression 


£g(r,  t)  - ^ 25  - ^ In^x  - j , 


a < r < b 


(19) 


Conservation  of  mass  in  the  elastic  region  of  the  first  layer  requires 

that 


o 


b < r < h 

o 


(20) 


The  approximate  relaticn  obtained  from  rearranging  Equation  20  and 
consolidating  the  resiilt  with  Equations  15  and  17  is 


£g(r,  t)  s 


b < r < h 
o 


(21) 
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(22) 


tic  state  ^ and  the  eir  confer  eat  ial  strain,  is 


jximated  "by 


EgCr,  t) 


3 

^2^0 


, "b  < h < r 
o 


which  is  analogous  to  Equation  21  for  the  first  layer.  The  subscript  "2® 
denotes  second  layer  material  properties.  Although  Equation  22  is  only 
approximately  correct ^ the  error  which  is  introduced  in  Equation  7 is 
the  same  order  of  mgnitude  as  the  error  which  results  from  the  assump- 
tions (a)  that  h ” h when  h < h and  (h)  that  the  entire  second 

o o 


layer  is  in  the  elastic  state  when  b < h 


Then  the  substitution  of 


b and  b^  for  r and  r^  , respectively 9 in  Equation  I8  results  in 


b < h 


(23) 


where 


6.  = a + 33. 


1*  2 


(2h) 


and  the  subscripts  ”l"  and  ”2”  again  denote  first  and  second  layer  mate- 
rial properties,  respectively,  Note  that  the  criterion  b < h is  now 

1/3  ° 

equivalent  to  a < 6^  ° 

11,  IJow  considering  the  situation  when  b > h li.e.s  when 

1/3  \ ® \ 

a ' ^oi®  Equations  18  and  19  are  applicable  to  the  locked  plastic 

regions  in  both  layers . Furthermore,  the  motion  of  the  interface  is  no 

longer  neglible,  and 


2 " 2 ® 

hh-aa  ,b>h 


Integrating  Equation  25  from  to  t yields 

a ^ (1  « 6. )h^  9 b > h 
1 o “ c 
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where  the  initial  cosdJ.tions  are  h(t  ) ~ h and  a(t  ) = 6, 

O O O pj. 


1/3, 


o 


le  cireumfereatial  strain  in  the  elastic  region  is  obtained  "by  replac- 
ig  6.  with  in  Equation  21  so  that 


Spb-" 

t)  a-»~  ^ 


r > b > h 


o 


(27) 


When  a"^  >>  , the  relation  between  a and  b must  approach 

3 3 1 ® 

b — a /6  , which  is  analogous  to  Equation  23.  Thus®  the  correspondini 

& 

d 9 

rate  equation  relating  a and  b must  approach 


^2 


(28) 


where  is  given  by  Equation  2Uo  In  order  to  obtain  an  approximate 

but  continuous  expression  for  b , Equation  28  is  integrated  from  t^ 
to  t j resulting  in 


■5  .3  “ 6 

^3  ^ ^ b > h 

^2  ° ° 


(29) 


1/3, 


where  the  initial  conditions  are  b(t  ) = h and  a(t  ) = 6.  "^h  . 

o O 0.0 


Mormal  Stress  at  the  Cavity  Surface 


12p  Approximations  have  now  been  established  relating  ^0(^s  "t)® 
b(t)  s h(t)  8 and  a(t)  for  the  two  cases  b < h and  b ^ h (i.e.j 
a < 6^  and  a 2.  “q  » respectively).  For  the  first  case. 

Equations  9?  11 ® l^s  19*  21,  22,  and  23  are  used  with  h — to 
evaluate  Equation  7.  For  the  second  case.  Equations  9»  11 » l^s  19?  26, 
27,  and  29  are  used  to  evaluate  Equation  7.  In  each  case,  the  resulting 
expression  for  the  compressive  normal  stress  at  the  cavity  surface  has 
the  form 


P = Pg  Pj  (30) 

where  p is  the  contribution  due  to  material  shear,  and  p^  is  the 

s * 


e 


front  at  the  initial  layer  interface  location ^ the 


egressions  for  p ^ B ^ and  are 


p (h  < h ) = ~ -I-  Y,  In  + “ E 
^ s o 3 1 1 9 ti 


/ o.  \ 


• I -a(‘  *a  (^' 


B^(b  < h^)  = 1 - 4/3  ^ 


02  \ a 


^01/  ^o 


BgCh  < h^) 


(1  - \r 

2B  -s-  ”2/0  — 

^ 6^'^(1  » a) 


I ^ 


V “1/  M3 

iirr-jT  ri 


^02  \ a"* 

"01/ 7 


(37) 


where  h and  b are  given,  respectively,  by  Equations  26  and  29  for 
b'  > h 

o 

13»  When  both  layers  have  identical  properties  (i.e.,  when  the 
medium  is  homogeneous),  the  layer  subscripts  are  dropped,  and  Equa- 
tions 31-37  reduce  to 

Pi  = 


(39) 


l6 


= 1 - 


g 

773 


1 ~ a J 2 


L - g/6)^ 
1 ” a 


(Ui) 


Material  Compressibility 


lU,  Throughout  the  development  of  the  foregoing  cavity  expansion 
theory  the  compaction  coefficient  a — - <<  1 has  been  regarded  as  a 

material  constant.  This  restriction  is  relaxed  by  allowing  a to  be  a 
slowly  varying  f\inction  of  pressure  and  then  neglecting  its  derivatives. 
By  taking  this  approach,  it  is  possible  to  relate  the  appropriate  in- 
stantaneous value  of  a to  the  dynamic  pressure  in  the  material  adja- 
cent to  the  cavity.  Since  the  theory  is  already  limited  to  small  values 
of  a , materials  of  interest  are  incompressible  to  a first  approxima- 
tion. For  an  incompressible  homogeneous  material,  the  volume-averaged 
dynamic  pressure^  in  the  plastic  region  a < r < b is  approximately 


P ^ P 
ave  2 


0 


6^3  . , , M3 

I I I > 


26j'^^  - 36j  + 


ik2) 


where 


= 3S  « 1 


(1^3) 


By  using  Equation  h2  in  conjimction  with  the  pressure-density  relation 
for  a given  material,  a first  approximation  can  be  obtained  for  the  com- 
paction coefficient.  When  is  sufficiently  small,  many  materials 

can  be  characterized  by  the  linear  pressure-density  relation 


«<  1 


(hh) 


where  C is  the  dilatational  wave  velocity  in  the  material. 

15.  The  cavity  expansion  theory  for  a concentrically  layered 


IT 


mediuai  is  restricted  to  materials  with  identical  compaction  coefficients. 
However,  this  theory  can  be  applied  to  materials  with  different  compac- 
tios  coefficients  by  estimating  a composite  compaction  coefficient  for 
both  layers  according  to 


a = 


where  and  obtained  from  Equations  U2  and  by  using  first 
and  second  layer  properties,  respectively.  The  effective  locked  plastic 
densities  for  the  two  layers  are  then 


P 


£1 


ih6) 


^£2 


(i?7) 


The  particular  expression  for  a which  appears  in  Equation  ^5  is 
chosen  so  that  the  composite  compaction  coefficient  will  vary  smoothly 
in  a manner  similar  to  the  inertial  coefficients  emd  . 
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PMT  III:  PHIOTATIOl  mOl' 


Baekgroiand 


2 3 

analogy  with  dynaaiic  spherical  cavity  expansion j Goodier  later  de= 

veloped  a theory  for  the  high-speed  penetration  of  incompressible, 

strain-hardening  targets  hy  rigid  spherical  projectiles.  Subsequently, 

Hanagud  and  Ross  modified  the  dynamic  spherical  cavity  expansion  theory 

in  an  effort  to  aecotmt  for  material  compressibility  and  incorporated 

the  results  in  Goodier ®s  penetration  theory.  The  next  development  came 

5 

when  Bernard  and  Hanagud  modified  the  Goodier  theory  in  an  attempt  to 
account  for  the  effect  of  projectile  nose  shape,  Bernard  and  Hanagud 
also  proposed  a dimensionless  parameter  called  the  "solid  Reynolds 
number"  which  serves  as  an  order-of-iiagnitude  index  of  the  ratio  of 
dynamic  pressure  to  shear  stress  in  the  target  and  which  can  be  used  to 
some  extent  to  determine  the  relative  applicability  of  the  theory  in  a 
particular  penetration  problem.  Moreover,  it  was  found  that  the  range 
of  applicability  for  the  penetration  theory  in  Reference  5 is  restricted 
to  low  to  moderate  values  of  the  solid  Reynolds  number  (corresponding 
to  shallow  to  moderate  penetration  depths).  At  high  values  of  the  solid 
Reynolds  number  (i,e.,  at  very  high  velocities)  the  theory  overpredicts 
the  target  resistance,  resiilting  in  an  underprediction  of  final  penetra- 
tion depth.  Accordingly,  the  Bernard-Hanagud  version  of  the  penetration 
theory  will  be  referred  to  as  the  "shallow  penetration  theory,"  and  the 
objective  of  the  present  Invistigation  will  be  to  develop  a deep  penetra- 
tion theory  which  is  applicable  for  very  deep  penetration  as  well  as  for 
shallow  penetration  in  cohesive  targets. 


Shallow  Penetration  Theorj 


17,  The  basic  assumptions  in  the  shallow  penetration  theory  are 
as  follows ; 
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_a.  Tlie  projectile  is  completely  rigid. 

The  portion  of  the  projectile  frontal  surface  which  is 
embedded  in  the  target  is  in  complete  contact  with  the 
target. 

Tangential  stresses  on  the  projectile  frontal  surface 
are  negligible. 

The  compressive  noirmal  stress  p on  the  frontal  surface 
is  given  by 

P - Pg  Pj  (^S) 


where  p is  the  contribution  due  to  target  behavior  in 
shear  an!  is  the  dynamic  pressure  in  the  target  ad- 

jacent to  the  frontal  surface. 

The  shear  contribution  p is  independent  of  projectile 
geometry*  uniform  over  the  frontal  surface,  and  equal  to 
the  cavity  expansion  shear  term,  which  is  given  by  Equa- 
tion 39  in  Part  II. 

l8.  With  these  assumptions  in  mind,  only  the  dynamic  presstire 
term  p^  remains  to  be  determined  in  the  expression  for  p . By  draw- 
ing a loose  analogy  with  the  cavity  expansion  theory.  Equation  26  is 
replaced  with 


for  fully  embedded  axisymmetric  projectiles  (Figures  3 and  ) , 
where 


(h9) 


= locked  plastic  density  of  the  target 
^1*^2  ~ coefficients  given  by  Equations  and  Ul 

= projectile  radius 

V = particle  velocity  in  the  target  adj acent  to  the  frontal 

^ surface 

q = projectile  acceleration  (the  latter  three  terms  replace 
a , a , and  a , respectively) 

This  then  leaves  the  distribution  of  the  particle  velocity  v tc  be 

P 

specified  along  the  frontal  surface. 
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Figure  3.  Projectile  with  fully  embedded 
conical  nose 


Figure 


Projectile  with  fully  embedded 
ogive  nose 
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19.  The  constraints  on  particle  motion  along  the  frontal  surface 

ai*e  as  follows  t 


a.  The  velocity  of  a particle  at  the  nose  tip  ^2  = 0 must 
he  identical  with  the  projectile  velocity  q . 

h.  Material  cannot  cross  the  projectile-target  interface; 
thus  the  particle  velocity  compc  nent  v , which  is 
normal  to  the  frontal  surface , must  be  equal  to  the 
normal  component  of  the  projectile  velocity. 

c_.  The  particle  velocity  must  he  continuous  on  the  frontal 
surface. 

ji.  Finally  it  is  assumed  that  the  particle  velocity  com- 
ponent v^  s which  is  tangent  to  the  frontal  surface, 
is  zero  at  the  base  of  the  nose  z = L , where  z is 
axial  distance  from  the  nose  tip  and  L is  projectile 
nose  length. 

20.  For  the  conical  nose  shown  in  Figure  3,  is  given  by 

v^  = q sin  (50) 

where  is  the  cone  half-angle . The  variation  of  v^  is  assumed  to 

be 


= (1 


.1/2  • 

x;  q cos 


<{>  , 0 ^ X ^ 1 


(51) 


where  the  nondimens ional  cone  position  is  denoted  by 


X 


(52) 


The  resulting  expression  for  v 


is  then 


V 

P 


2 

X cos 


For  the  ogive  nose  shown  in  Figure  U, 


(53) 


V = q sin  n 
n 


(5M 


and  the  equivalent  dimensionless  cone  position  is 


5 
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(55) 


x(^) 

e 


cos  n 


cos  n - cos  4> 

1 - cos  4<  cos  n 


where  Equation  55  is  analogous  to  Equation  52.  Replacing  x and  <|) 
with  n » respectively,  in  Equation  51  leads  to 


1/2  • 

[1  --  x^(n)]  q cos  n » 0 _<  n _<  (}> 


(56) 


so  that 


V 

P 


[1  - X (n)  cos^ 

e 


(57) 


21.  To  obtain  the  projectile  equation  of  motion  for  the  case  f 
a conical  nose.  Equation  53  is  inserted  into  Equation  ^^9*  and  the  re- 
sulting expression  for  is  used  in  Equation  48.  The  total  axial 

resisting  force  F acting  on  the  projectile  is  then  fovind  according 

to 


F 

z 


[Pg  + Pj(Z, 


<{))  ] Sin  (f)  dA(z,  (j>) 


(58) 


where  the  surface  element  dA(z,  4>)  is  expressed  as 

dA(z,  = 2tt  z dz  (59) 

cos  9 

Similarly,  tor  an  ogive  nose.  Equation  57  is  inserted  in  Equation  49, 

and  the  resulting  expression  for  is  used  in  Equation  48.  Then 

F is  found  according  to 
z 


I' 

z 


[p^  + Pt(h)j  sin  n dA(n) 

3 X 


(60) 


where  the  surface  element  dA(n)  is  expressed  as 
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(61) 


dA(n) 


2iTa  (cos  n - cos 

o 


cos  ^ - cos  r\)  cos  ti  - cos  r|]  dr) 


(l  -■  cos  4>)  cos  n 


For  both  cases,  the  projectile  equation  of  motion  is  obtained  by  setting 

• • 

Mq  = ~F  , where  M is  the  projectile  mass.  The  resulting  equation  of 
motion  has  the  form 


where  p , B , and  B are  given  by  Equations  39-^1  as  functions  of 
the  target  yield  strength  Y , elastic  modulus  E , strain-hardening 
modulus  , and  compaction  coefficient  a . The  quantity  f^  is  the 
so-called  nose  shape  factor  resulting  from  the  variation  of  dynamic 
pressure  along  the  projectile  frontal  surface.  For  conical  noses,  f^ 
is  given  by 

f = sin^  4*  + T cos^  <j)  (63) 

n 3 


where  ^ is  related  to  the  nose  length  L and  projectile  diameter 

D = 2a  by 
o 


tan  4>  = ^ 


(6M 


For  ogives,  the  nose  shape  factor  is 

f = 1 „ - ^ B^  In  (2  - e)  - (3B^  + 2B)(B  - e) 

^ - e)^  I 

+ j (3B^  + 6B  + 1)(B^  - e^)  _ 1 (B^  _ 

- j (B^  -s-  6b  + 3)(B^  - ^ 1 ^ 


(65) 


2k 


where 


3 = 2e  - e 


(66) 


and  the  parameter  e 
and  to  L/D  by 


is  related  to  the  ogive  caliber  radius  (CRH) 


CRH  = 


R 

D 


2^ 

2e 


(67) 


(A  hemisphere  is  an  ogive  with  CRH  = L/D  = 1/2  , in  which  case  Equa- 
tion 65  reduces  to  = 2/3  . ) As  pointed  out  in  Reference  5»  Equa- 
tion 62  is  pre'^’omed  to  be  valid  only  for  projectiles  and  targets  for 

3 

which  the  condition  ira  <<  M is  satisfied.  In  any  case,  the  term 

o S,  1 -V  j 

is  thought  to  represent  an  upper  limit  insofar  as  the  effect  of 
target  acceleration  is  concerned. 

22.  Daring  the  embedding  of  the  nose  (penetration  depth  q £ L), 

the  projectile  equation  of  motion  is  obtained  by  substituting  a for 

a in  Equation  62,  where  a is  the  cylindrical  radius  of  the  pro- 
o s 

Jectile  ^ the  target  surface. 

23.  For  the  penetration  of  layered  targets,  the  quantities  p , 


1 


, and  Bg  are  modified  in  accordance  with  the  cavity  expansion  theory 


for  a concentrically  layered  medium.  The  derivation  of  these  quantities 
in  Reference  5 is  lengthy  and  cumbersome;  thus  it  is  appropriate  that 
the  discussion  for  layered  targets  be  deferred  \ontil  later  when  a 
simpler  approach  will  be  taken  in  the  development  of  the  deep  penetra- 
tion theory. 

Solid  Reynolds  Number 

2k.  Among  the  various  parameters  defined  in  Reference  5 is  the 
"solid  Reynolds  mmber,"  which  is  an  or der-of -magnitude  index  of  the 
ratio  of  dynamic  pressiize  (i.e. , inertial  stress)  to  shear  stress  in 
the  target,  defined  by 


R 

^s  ~ Y 


:68) 
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The  solid  Reyiiolds  number  can  be  used  as  a scaling  parameter  for  ex- 
perimental and  theoretical  results  for  cohesive  targets  and  also  es  an 
indicator  of  the  relative  applicability  of  the  shallow  penetration 
theory  to  particular  penetration  problems.  The  working  range  of  this 

theory  for  depth  predictions  appears  to  be  0 < R 5 100  ^ while  the 

s 

working  range  for  projectile  deceleration  predictions  (accurate  within 

a factor  of  two ) is  more  like  0 < R i 10  . In  other  words , the  shal- 

s 

low  penetration  theory  seems  to  work  best  for  cases  in  which  the  re- 
sistance of  the  target  due  to  shear  is  comparable  to  that  due  to 
inertia.  When  the  solid  Reynolds  number  becomes  very  large  (i.e. , for 
"high"  velocities  and  "soft"  targets ) , the  theory  overpredicts  the  pro- 
jectile deceleration,  which  results  in  an  underpredicticn  of  the  final 
penetration  depth.  Furthermore,  the  magnitudes  of  the  overprediction 


and  underprediction  tend  to  increase  with  increasing  R^  . 

25.  The  breakdown  of  the  shallow  penetration  theory  at  high 

values  of  R is  not  siirprising  since  the  dynamic  pressxire  variation 
s 

along  the  projectile  frontal  surface  is  obtained  from  a particle  veloc- 
ity  distribution  v /q  , which  is  a fvmction  only  of  projectile  geometry 
and  axial  distance  from  the  nose  tip.  This  resiilts  in  a nose  shape 
factor  or  "drag  coefficient"  which  is  a function  of  projectile  geometry 
alone  and  which  is  independent  of  projectile  velocity  and  target  proper- 


ties. In  order  to  account  for  the  variation  of  with  projectile 

geometry,  projectile  velocity,  and  target  properties , an  approach  is  now 
sought  that  will  yield  as  a fiinction  of  the  solid  Reynolds  number 

as  well  as  the  projectile  geometry. 


Penetration  Theory  for  Homogeneous  Targets 


26.  In  attempting  to  develop  a penetration  theory  which  is  appli- 
cable at  high  values  of  R , one  might  modify  any  or  all  of  the  basic 

s 

assumptions  of  the  shallow  penetration  theory  (see  paragraph  17 ) • In 
particular,  the  assumption  of  complete  contact  between  the  target  and 
the  projectile  frontal  surface  is  questionable , and  the  possible  occur- 
rence of  separation  on  the  frontal  surface  may  affect  the  resistance  of 
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the  target  to  penetrationo  Nevertheless,  all  of  the  basic  assumptions 
of  the  shsillow  penetration  theory  vill  "be  retained,  and.  the  development 
of  the  deep  penetration  theory  will  proceed  from  a simple  modification 
of  the  particle  velocity  distribution  on  the  projectile  frontal  surface » 
2T»  The  constraints  on  particle  motion  along  the  frontal  surface 
(see  paragraph  19)  are  taken  to  be  the  same  as  in  the  shallow  penetra- 
tion theory,  except  for  the  assumption  that  v^  ~ 0 at  the  base  of  the 
nose.  This  particular  constraint  is  relaxed  to  allow  more  flexibility 

in  the  distribution  of  v.  along  the  frontal  surface, 

t 

28„  Considering  first  a projectile  with  a fully  embedded  conical 
nose  as  shown  in  Figure  3,  the  particle  velocity  component  « normal 
to  the  frontal  stirface,  is  still  given  by  Equation  50,  However,  the 
tangential  velocity  component  v^  is  now  assumed  to  vary  according  to 

-4»(R3)x^/2 

V (x,  R ) = (q  cos  <i>)e  (69) 

u S 


« 2 

where  x = z/L  , R = p_q  /Y  , and  t|;(R  ) is  an  unknown  function  of 

s u s 

R . Equation  69  is  by  no  means  a unique  expression  for  v . It  does, 

S u 

however,  represent  a convenient  fimctional  form  in  which  v^  decreases 
with  increasing  x and  increasing  i|<  . If  <|'(R  ) is  chosen  to  be  a 

5 

monotonically  increasing  function,  then  Equation  69  will  result  in  an 
expression  for  the  nose  shape  factor  or  "drag  coefficient"  that  de- 
creases with  increasing  projectile  velocity  and  nose  length  (which  is 
more  or  less  consistent  with  experimental  observation).  Incorporating 
Equation  69  (which  now  replaces  Equation  51 ) in  Equation  53,  the  expres- 
sion for  V becomes 
P 


V 

P 


sin 


■¥  cos 


2 


(70) 


Equation  70  is  substituted  into  Equation  which  is  in  turn  inserted 
into  Equation  58 » Setting  Mq  = -F  as  before,  the  resulting  projec- 

Zi 

tile  equation  of  motion  is 
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where 


4-  na"  p„l 
o i 


l)«  = -”“o 


f^(d).  R^)  = sin^  4*  + cos^  ^ (T2) 

s 

Equation  12  now  replaces  Equation  63  for  the  nose  shape  factor 

during  all  phases  of  the  penetration  process,  including  the  embedding  of 

the  nose, 

29.  Equation  12  is,  of  course,  restricted  to  projectiles  with 
conical  nose  shapes.  In  order  to  apply  the  present  approach  directly 
to  ogives.  Equation  69  should  be  replaced  by 

-'{'(R  )x^(n)/2 

V.  (n,  R_)  = (q  cos  n)e  (73) 

T>  S 


where  n is  defined  in  Figure  ^4,  and  x (n)  is  given  by  Equation  55. 

® 2 2 2 

Equation  73  should  then  be  used  in  the  relation  v * v + v,  which  is 

p n t 

in  turn  substituted  into  Equation  k9.  The  resulting  equation  for  the 
dynamic  pressure  p^  must  then  be  inserted  into  Equation  60  to  obtain 
the  resisting  force  F and,  subsequently,  the  projectile  equation  of 
motion  which  will  be  of  the  same  form  as  Equation  71  but  with  an  ex- 
tremely complicated  relation  among  f^  , , end  the  maximum  half- 

angle 4>  (Figure  4).  Moreover , the  integral  for  F must  be  evaluated 
numerically  when  this  approach  is  taken,  and  thus  it  is  expedient  to 
seek  a simpler  approach  that  results  in  an  approximate  analytic  expres- 
sion for  f 


n 

30.  In  order  to  obtain  an  expression  for  f^  which  has  the 
SMie  form  as  Equation  72,  the  effective  ogive  half-angle  is  defined 
as  follows : » 
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The  nose  shape  factor  for  ogives  is  then  approximated  by 

f (R  , ♦ ) ^ sin^  ♦ + "~-77^1 cos®  4'  (78) 

n s ^e  ^e  ij;(R  ) e 

s 

According  to  this  scheme,  any  ogive  can  be  replaced  by  its  equivalent 
cone  which,  by  definition,  has  approximately  the  same  value  of  f^  as 
the  ogive o Equivalent  cone  L/D  ratios  are  listed  for  ogives  in 
Table  1. 

31.  At  this  point  only  the  relation  between  and  R is  still 

s 

unspecified.  In  accordance  with  paragraph  28,  ip  is  defined  as  a mono- 

tonically  increasing  function  of  R with  the  limiting  values 

s 


29 


(79) 


11m  ip  = 0 

R 0 


lim  4;  - «} 

R oo 

s 


(80) 


Table  1 

Equivalent  Cone  Parameters  for  Ogives 


Ogive  CRH 

Ogive  L/D 

Equivalent 
Cone  L/D 

Ogive  CRH 

Ogive  L/D 

Equivalent 
Cone  L/D 

0.5 

0.50 

0,50 

5.5 

2.29 

1.99 

1,0 

0.8t 

o.t8  : 

6.0 

2.U0 

2.08 

1.5 

1.12 

0.99 

6.5 

2.50 

2.17 

2.0 

1.32 

l.l6 

7.0 

2.6o 

2.26 

2.5 

3.50 

1.31 

7.5 

2.69 

2. 31+ 

3,0 

1.66 

1.1+5 

8.0 

2.78 

2.1+2 

3.5 

1.80 

1.57 

8.5 

2.87 

2.50 

H.O 

1.9i» 

1.69 

9.0 

2.96 

2.57 

2,06 

1.8o 

9.5 

3.0I+ 

2.61+ 

5.0 

2.18 

1,90 

10.0 

3.12 

2.71 

The  corresponding  limiting  values  for  are  then 

lim  f =1 
n 

R -»■  0 
s 


(81) 


with 


2 

lim  f - sin  4* 
n 

R “ 
s 


(82) 
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for  cones  and 


2 

lim  f - sin 
n 

R -4-00 

s 


(83) 


for  ogiYes.  It  now  remains  to  determine  a single  monotonically  increas- 
ing function  i^(R  ) which  satisfies  Equations  T9  atid  80  and  which  re- 

s 

suits  in  acceptable  predictions  of  projectile  deceleration  and  final 
penetration  depth,  TOiis  function  must  be  determined  empirically ^ but 

once  established  it  will  be  considered  "universal,”  i ,e. , it  will  be 

presumed  to  be  the  same  for  all  projectiles  and  all  cohesive  targets, 

within  the  context  of  the  present  theory.  If  a different  function  ) 

s 

has  to  be  chosen  for  each  target  and  each  projectile,  then  the  theory 

is  no  better  than  a strictly  empirical  approach  in  which  the  "constants" 

in  the  projectile  equation  of  motion  must  be  determined  separately  for 

each  individual  situation.  On  the  other  hand,  if  a single  function 

if/(R  ) can  be  found  that  correlates  experimental  data  for  a variety  of 
s 

projectiles  and  targets , for  both  shallow  and  deep  penetration,  then 
the  theory  may  represent  a usefxil  tool  for  making  predictions  before  the 
fact, 

32.  With  the  preceding  observations  in  mind,  the  following  equa- 
tion is  proposed  for  i|;(R  ) : 

s 


It  will  be  shown  in  Part  IV  that  Equation  results  in  good  agreement 

between  theory  and  experiment  for  a number  of  different  projectiles  and 

targets.  Other  expressions  for  if>(R  ) have  been  found  that  also  produce 

s 

good  results  (e.  g. , = R /15 ) , but  Equation  8U  appears  to  produce  the 

s 

best  overall  results  within  the  context  of  the  present  theory. 


Deep  Penetration  Theory  for  Layered  Targets 
33.  In  order  to  modify  the  projectile  equation  of  motion  (71 ) for 
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the  penetration  of  layered  targets,  the  appropriate  expressions  for  the 
quantities  p , B , and  B must  first  he  determined  hy  means  of  an 

S J.  ci. 

analog  with  dynamic  cavity  e::^ansion  in  a concentrically  layered 
medium.  The  situation  to  he  considered  is  illustrated  in  Figure  5; 


Figure  5.  Projectile  penetration  in  a layered  target 


however,  the  target  may  have  any  number  of  layers  as  long  as  the  pro- 
jectiie  motion  is  significantly  influenced  by  only  two  layers  at  a time. ' 
The  penetration  process  is  now  divided  into  two  regimes:  (a)  when 

L < H , the  nose  tip  has  not  yet  reached  the  layer  interface;  and 
(b)  when  L ^ H,  the  projectile  has  begun  to  perforate  the  layer  inter- 
face. Each  of  these  regimes  can  be  divided  into  three  subregimes: 

(a)  when  S < H , the  plastic  front  has  not  yet  reached  the  layer 
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terface;  and  (c)  when  H , the  effect  of  the  first  layer  on.  the 

projectile  equation  of  aotion  is  negligible® 

3h.  To  approximate  the  values  of  Pg  s ® 
above  regimes  and  suhregimeSs  it  is  first  necessary  to  consider  Equa- 
tions 30-^1  in  Part  II®  In  Reference  5s  an  attempt  is  made  to  account 

« 

for  the  deformation  of  the  layer  interface  due  to  the  approaching 

projectile.  Moreover , en  upper  limit  approximation  is  used  to  obtain 
the  effective  interface  deformation,  and  the  i-esulting  expressions  for 
H , S , p g g and  B_  are  not  only  lengthy  and  cumbersome,  but 
the  approximation  for  the  interface  deformation  is  itself  questionable 
and  has  only  a small  effect  on  the  actual  values  of  Pg  s s 
A great  deal  of  complication  can  then  be  avoided  (without  any  demon- 
strable loss  of  acctsracy)  by  simply  neglecting  the  deformation  of  the 
interface  altogether . Accordingly,  the  quantities  h(t)  and  h which 
appear  in  Equations  32-37  will  be  replaced  by  the  distance  H from 
the  base  of  the  projectile  nose  to  the  interface, 

H = H - q -f  L (85) 

o 

As  a result,  the  criteria  L < H and  L ^ H can  now  be  replaced  by 

q < and  q ^ » respectively. 

35.  In  the  first  regime  where  q < H , the  equations  for  p , 

o s 

B,  5 and  B_  are  obtained  by  replacing  a(t)  with  a , h(t)  and  h 
Id  o o 

with  H 3 and  b(t)  with  S (effective  distance  from  base  of  projec- 
tile nose  to  plastic  front ) in  Equations  23,  29,  and  32-37*  Then , 

S < H is  substituted  for  the  criterion  b < h . The  restilts  are  as 

o 

follows : 

^3 

^ , S < H (86) 
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Thus,  when  q < H , p is  substituted  for  p in  the  projectile 

o ^ 

equation  of  motion  (7l),  and  the  expressions  for  p , B , and  B are 

S J.  ^ 

replaced,  respectively,  by  Equations  87-89  when  S < H and  by  Equa- 
tions 91-93  when  S _>  H . The  quantity  a is  the  composite  compaction 
coefficient  for  the  two  layers  (which  will  be  discussed  later);  6^  and 
6^  s-re  defined  by  Equation  2U;  and,  with  the  exception  of  B^  and  B^  , 
the  subscripts  "l"  and  "2"  denote  quantities  evaluated  in  the  first 
and  second  layers,  respectively.  In  the  regime  q < , the  effective 

nose  shape  factor  f^.j^  is  approximated  according  to 


(9M 


. 2 

where  f , and  f ^ are  obtained  by  inserting  R ^ q /Y^  and 

nig  n2  si  01  ^ 1 

^s2  ~ *^02^  ’ respectively,  into  the  equation  for  f^  . Thus,  when 

q < H , then  f = f * and  the  product  f B^  is  given  by 
^o  nn*  ^ n2° 


f B„  = f *B_(S  < H)  , S < H 
n 2 n’'  2 


(95.) 


^n®2  = ^n^®2^^  - , S > H 


(96) 


When  a = H , Equation  9^  reduces  to  f * = f _ , Equation  90  becomes 
o n’^  n2 
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approximately 


3 

q ^ 

S ^ » a ^ H 
5 ^ * o 


(97) 


and  Equations  8T"89  and  91-93  give  way  to  Equations  39-^1 ■ (That  isj 
the  projectile  equation  of  motion  is  the  same  as  for  a homogeneous 
target  which  has  the  same  properties  as  the  second  layer )»  According?-y . 
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Pp  is  replaced  in  Equation  71  hy  and  p 
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36.  In  the  other  regime  where  q H , the  expressions  for  p ^ 

' o s 

. and  f when  S < H are  approximated  by 
1 n 2 
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where  is  the  cylindrical  radius  of  the  projectile  ^ the  interface 

location  E ^ and  S is  still  given  "by  Equations  90  and  95^  The  con- 
ditions  S < H and  q > H can  occur  simultaneously  only  when  L > S . 
i.e.j  'When  the  projectile  nose  length  is  greater  than  the  effective 
radius  of  the  plastic  front  relative  to  the  base  of  the  nose.  On  the 
other  hand,  when  S ^ H , » ®1  * ^n®2  approximated  by 


P 


s 


a 


p^(S  > H)  + -|  p^(K  < a^) 

a 

c 


P AQ^T 
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(106) 


For  both  regimes  (i.e.,  q < and  9L  ^ "the  expi'essions  for  , 

, and  reduce  to  Equations  98-100  when  = H . Thus,  whenever 

H j<  , the  projectile  equation  of  motion  is  Equation  71  with 
Po“Poo*  ^ P » ®P  given  by  Equations  98- 

100,  respectively. 


Target  Compressibility 

37.  In  order  to  estimate  the  target  compaction  coefficient  a 
at  any  time  during  the  penetration  process,  an  estimate  is  first  made 
for  the  volume-a /eraged  dynamic  pressure  plastic  zone  sur- 

rounding the  projectile  frontal  surface.^  This  is  done  by  substituting 
q for  a in  Equation  h2  and  dropping  the  acceleration  term  (projectile 
acceleration  is  always  negative  and  could  result  in  a negative  value  of 
^ave  * inadmissible).  The  resulting  equation  is  multiplied  by 

f to  account  for  the  effect  of  nose  shape  at  high  values  of  , and 
n s 

one  then  obtains 
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(107) 
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This  equation  is  used  in  a pressure~density  relation  (eog^^  Equation  hk] 
to  estimate  a ^ For  layered  targets,  a is  calculated  separately  for 
each  layer,  and  then  a composite  value  is  determined  "by  substituting  a 

o 

for  a(t)  and  H for  h(t)  in  Equation  ^5  yielding 

Q>  \ 3* 

a - * ^2  ^ 9 % 1 ® (108) 


38,  lo  modifications  have  been  made  in  the  oblique  impact  theory, 
which  is  discussed  in  Reference  The  oblique  theory  is  still  limited 
to  predictions  of  projectile  rotation  during  the  embedding  of  the  nose 
and  is  yet  to  be  verified  by  comparison  with  experiment.  However,  in 
light  of  the  present  work,  it  is  appropriate  to  replace  the  value  of  f 

n 

used  in  Reference  5 with  that  given  by  Equation  72.  Thus,  the  equations 
which  govern  the  forces  and  moments  acting  on  the  projectile  will  have 
the  same  form  as  before,  but  the  incorporation  of  Equation  72  for  f^ 
introduces  an  implicit  dependence  upon  the  solid  Reynolds  number.  Since 
f is  a decreasing  function  of  R , maximum  projectile  rotation  will 

II  s 

be  seen  at  low  values  of  R^  , and  the  projectile  rotation  will  be  a 

more  rapidly  decreasing  function  of  impact  velocity  than  that  produced 

by  the  value  of  f which  was  used  in  Reference  5. 

n 

Hondimensional  Results  for  Typical  Targets 


39-^  Although  the  penetration  theory  which  has  been  developed  may 
appear  somewhat  complicated  at  first,  the  implications  for  homogeneous 
targets  are  fairly  simple  when  viewed  in  the  proper  perspective.  For 
heavy  projectiles  (i.e.,  ira^p^B^  <<  M) , the  projectile  equation  of 
motion  reduces  approximately  to 
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For  slightly  compressible  targets s the  range  of  values  for  P£®2 
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and  the  corresponding  range  for  p is 
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It  is  reasonable  to  characterize  a "typical”  target  by  p = kY  and 
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p B q /p  = — R , so  that  Equation  110  becomes 
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When  the  transformation 
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are  introduced.  Equation  113  gives  vay  to 


Tta  p 

os 


dR 

dT 


dR 


VSTTa  f 


* 3 ‘■n«s  > 1 ^ 


(ll6) 


^ Sl+:|-fR  , q>L 
dC  3 n s * ^ — 


For  deep  penetration  (final  penetration  depth  q^  >>  L),  Equation  ll6 
can  be  replaced  by 
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The  peak  nondimen sional  deceleration  of  the  projectile  can  thus  be  ob~ 

2 

tained  by  using  the  impact  value  of  R = p^v  /Y  , and  the  corresponding 

5 U O 2 

dimensional  value  is  obtained  by  multiplying  by  Uira^Y/M  . The  final 
nondimensional  penetration  depth  = 8'>ra^PQqj./M  is  found  by  inte- 
grating Equation  117  niomerically,  and  the  corresponding  dimensional 
depth  is  found  by  multiplying  C^.  by  M/Sttb^Pq  . In  Figure  6,  peak 
nondimensional  deceleration  is  plotted  against  nondimensional  impact 
velocity  v^^’p^Ty”  for  conical  noses  with  L/D  = 1,  2,  and  3 . In 
Figure  7,  nondimen s i onal  final  penetration  depth  is  plotted  against  non- 
dimensional  impact  velocity.  Results  are  presented  for  the  shallow 
penetration  theory , in  which  = f^(<j>)  , and  for  the  deep  penetration 

theory,  in  which  f = f (4»,  R ) . For  very  blunt  noses,  i.e. , 

n n s 

L/D  ^ 1/2  , the  predictions  of  the  deep  penetration  theory  reduce  essen- 
tially to  those  of  the  shallow  penetration  theory.  Thus,  it  is  empha- 
sized that  at  high  values  of  R^  deep  penetration  theory  is  probably 
valid  only  for  sharp  projectiles  where  L/D  ^ 1 . A deep  penetration 


Figure  6.  Nondimen sional  peak  deceleration 
versus  impact  velocity  for  conical  projec- 
tiles and  ”t3rpical"  targets 


Figure  ?•  Wondimensional  final  penetration 
depth  versus  impact  velocity  for  conical 
projectiles  and  "typical”  targets 


theory  for  'blunt  projectiles  should  account  for  the  effect  of  target 
material  just  ahead  of  the  projectile  which  is  formed  into  a sharp  coni 
cal  "nose”  during  the  impact  phase. 


k2 


PART  IV:  APPLICATIOIS 


Introduction 


^*0,  This  part  of  the  report  is  devoted  to  the  application  of  the 
deep  penetration  theory  to  problems  of  practical  interest  and  to  the 
comparison  of  theoretical  results  with  experimental ^ empirical,  and  2D 
finite-difference  results  for  cohesive  targets.  In  the  applications 
which  follow,  the  range  of  the  solid  Reynolds  number  extends  from  very 
low  values  (shallow  penetration)  to  very  high  values  (deep  penetration), 
and  the  function  's’(R  ) is  given  by  Equation  8U  for  all  targets  and  all 
projectiles.  In  each  example,  the  predictions  of  the  shallow  penetra- 
tion theory  are  compared  with  those  of  the  deep  penetration  theory. 

Wien  comparisons  are  shown  with  2D  finite-difference  calculations,  no 
attempt  is  made  to  assess  the  validity  or  accuracj'  of  the  various  as- 
sumptions and  techniques  employed  by  different  authors  since  this  lies 
beyond  the  scope  of  the  present  investigation.  Thus,  the  reasons  for 
agreement  or  disagreement  between  theoretical  end  2D  solutions  remain 
open  to  speculation.  The  agreement  achieved  between  theoretical  predic- 
tions and  empirical  and  experimental  resizlts  is  encouraging  but  not 
conclusxve.  Furthermore,  the  relative  applicability  or  inapplicability 
of  the  deep  penetration  theory  to  the  overall  range  of  practical  problems 
can  only  be  demonstrated  by  extensive  comparison  with  experimental  data 
from  deep  penetration  tests.  However,  most  data  which  are  now  available 
lie  in  the  range  of  shallow  to  moderate  depths  for  which  the  theory  is 
already  fairly  well  doctunented. 

Choice  of  Material  Properties 

i+1.  The  development  of  the  deep  penetration  theory  stems  from 
the  cavity  expansion  theory  which  is  based  upon  the  presiamption  of  rate- 
independent,  bilinear  material  behavior  in  shear.  Unfortunately,  real 
juaterials  exhibit  neither  bilinear  shear  behavior  nor  rate- independence, 
and  so  there  exists  the  problem  of  choosing  an  appropriate  bilinear 


SHEAR  STRESS 


stress-strain  curve  that  adequately  approximates  real  material  “behavior 
insofar  as  it  relates  to  cavity  expansion  and  projectile  penetration. 

For  upper-“botind  predictions  of  final  penetration  depth,  the  results  of 
miconfined  static  trir.'^ial  compression  tests  are  used  to  estimate  target 
shear  properties;  and  for  many  applications,  these  upper-“bound  predic- 
tions are  adequate.  However,  for  ^ priori  predictions  in  untried  tar- 
gets, upper-hound  predictions  should  “be  supplemented  with  lower-hound 
predictions  using  properties  obtained  from  the  results  of  dynamic  tests. 
Anisotropic  behavior  may  further  complicate  the  situation  so  that  upper- 
and  lower-bound  properties  have  to  be  chosen  with  regard  to  maximum 
variation  due  to  anisotropy  as  well  as  rate  effects.  In  any  case,  ac- 
ceptable results  are  usually  obtained  by  setting  = 0 and  idealizing 
the  "real"  stress-strain  cxirve,  as  shown  in  Figure  8. 


Figure  8.  Idealization  of  a 
’’real"  stress-strain  curve 


Shallow  Penetration  of  Rock 


k2,  Thigpen^  has  performed  2D  finite-difference  calculations  in 

an  attempt  at  simulating  two  rock  penetration  tests  reported  by 
7 

Patterson.  In  these  calculations  the  shallow  penetration  of  nonrigid 
frietionless  projectiles  in  Madera  limestone  and  welded  tuff  is  analyzed 


identical  projectiles.  For  the  range  of  dynamic  pressures  encountered ^ 
the  compressibility  of  "both  targets  is  adequately  characterized  by 
Equation  and  the  elastic  raodtilus  is  approximately  equal  to  the  ini- 
tial value.  In  the  limestone  calculation*  the  target  properties  used 
by  Thigpen  are 

Pq  = 168  pcf®  (2.69  ^/cm^) 

Y = 13*690  psi  (0.9^^  kbar) 

E — 3.15  ^ 10^  psi  (217  kbar) 

E , = 0 

t 

C - 11,155  fps  (3^00  m/sec) 
and  the  (rigid)  projectile  characteristics  are 

W = 6lh  lb  (306  kg) 

D = 8 in.  (20.32  cm) 

CRH  = 9-25  (ogive  nose*  L/D  = 3) 

V = 5T0  fps  (IT^  m/sec) 
o 


where 

Pq  = initial  density 
Y = compressive  yield  strength 
E = elastic  modulus 
E “ strain-hardening  modulus 

Xt 

C = dilatational  wave  velocity 
W = weight 
D = diameter 

CRH  “ ogive  caliber  radf-us 


* A table  of  factors  for  converting  U.  S.  customary  units  of  measure- 
ment to  metric  (Sl)  units  is  presented  on  page  3. 


In  the  welded  tuff  calculation,  the  target  properties  used  hy  Thigpen  are 

pQ  " 115  pcf  (1  = 85  ^/m^) 

Y = 5510  psi  (0380  kbar) 

E S 1.09  X 10^  psi  (75.5  khar) 

^ 0 
t 

C - 6725  fps  (2050  m/sec) 


and  the  (rigid)  projectile  characteristics  are 

W = 1000  lb  (U55  kg) 

D = 9 in.  (22.9  cm) 

CRH  - 6.0  (ogive  nose,  L/D  ^ 2. h 

V = 695  fps  (212  m/sec) 
o 


The  solid  Reynolds  numbers  at  impact  are  R =0.86  and  2.17  in  the 

S 

limestone  and  tuff,  respectively.  The  results  are  compared  in  Figures  9 

and  10,  showing  that,  for  low  values  of  R , the  deep  penetration  the- 

5 

ory  predicts  slightly  higher  decelerations  and  shallower  final  penetra- 
tion depths  than  the  shallow  penetration  theory,  (This  trend  is  reversed 
at  high  values  of  » ) As  pointed  out  in  Reference  6,  the  oscillations 
in  the  2D  results  are  due  to  projectile  nonrigidity  which  is  accounted 
for  in  neither  the  shallow  nor  deep  penetration  theories.  No  experi- 
mental deceleration  record  is  available  for  the  welded  tuff  test. 

Shallow  Penetration  of  a Hypothetical  Frozen  Soil 

, 8 

43.  Ito  et  al.  have  performed  a 2D  finite-difference  calculation 
analyzing  the  penetration  of  a hypothetical  frozen  soil  by  a rigid  pro- 
jectile. The  2D  calculation  incorporates  an  elastic-plastic  material 


DECELERAl 


D FINITE-DIFFERENCE  RESULT 


Ka  7 Fil  4 ^RliWai  4 k 1 4 il^  I 


DEEP  PENETRATION  THEORY 


rL4  ■■si  4^11  Ji 


OBSERVED  FINAL  DEPTH® £.®®  FT 


0 O.IS  0.30  0.4S  o.eo  0.75  < 

DEPTH, M 

1  I I I I I 

0 0.5  1.0  1.5  ^.0  2.5 

DEPTH,  FT 

Figure  9-  Projectile  deceleration  versus  depth 
in  Madera  limestone 


2D  FINITE-DIFFERENCE  RESULTS 

SHALLOW  PENETRATION  THEORY 

DEEP  PENETRATION  THEORY 

OBSERVED  FINAL  DEPTH=4.95  FT 


W=IOOO  L0 
0 = 9 IN. 

CRH  =6 

Vq=695  ft/sec 
I I 

o.e  1.2 

depth, m 

4 

DEPTH, FT 


Figure  10.  Projectile  deceleration  versus  depth 

in  welded  tuff 


The  unconfiiied  target  properties  are 

pp  “ 125  pcf  (2,00  ga/cm^) 

Y ™ 2030  psi  (0,1^0  kbar) 

E ^ 5. TO  « 10^  psi  (39.3  kbar) 

™ 0 
t 

and 

C — TlOO  fps  (2165  m/sec) 


Taxget  compressibility  is  adequately  characterized  by  Equation  UU  for 
the  range  of  dynamic  presstires  encountered.  The  projectile  character- 
istics are 

W = 1T.3  lb  (T086  kg) 

D = 3 in.  (T.62  cm) 


and 


CRH  = 2.35  (ogive  nosej  L/D  = I.U5) 
v^  = U50  fps  (13T  m/sec) 


The  shallou"  and  deep  penetration  theories  are  now  applied  to  the  same 

problem  (in  the  absence  of  friction) 9 and  the  results  are  compared  in 

Figures  11  and  12.  The  solid  Reynolds  number  at  impact  is  R = 2.69  9 

s 

and  both  theories  produce  peak  decelerations  and  final  penetration 
depths  that  are  in  approximate  agreement  with  the  2D  results  for  the 

"frictionless"  case. 

Shallow  Penetration  of  Alumintun 
44.  Since  the  eavity-e^ansion-based  penetration  theory  has 
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DECELERATION 


Figure  11.  Projectile  deceleration  versus  time 
in  hypothetical  frozen  soil  target 


Figure  12.  Penetration  depth  versus  time 
in  hypothetical  frozen  soil  target 
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changed  somewhat  since  Its  initial  formulation  "by  Goodier,  it  is  appro- 
priate to  apply  the  present  shallow  and  deep  penetration  theories  to 
at  least  one  of  the  problems  which  were  discussed  by  Goodier  in  his 
original  work*  Consider  the  penetration  of  a homogeneous  202i+-T3  alumi- 
target  by  sphei^ical  steel  projectiles.  The  target  properties  are 

pQ  = 173  pcf  (2.77  gm/cra^) 

Y = 5.6  X 10^  psi  (3.86  kbar) 

E = lo"^  psi  (689  kbar) 

= li.l  X 10^  psi  (28.3  kbar) 

and 

— 10^  psi  (689  kbar) 


Target  compressibility  is  adequately  characterized  by  Equation  for 
the  range  of  dynamic  pressures  encountered.  Goodier  gives  his  results 
in  nondimensional  form,  choosing  2q^/D  as  the  nondimens, ional  final 


penetration  depth  and  v ^p  /Y  as  the  nondimensional  impact  velocity. 

o p 

The  quantity  p represents  the  projectile  density,  where  p = ii91  pfiY 
■D  P P 

(7.82  gm/cm^)  for  steel.  Consequently,  a nondimensional  impact  velocitjr 

v^'^Ij^Ty”  = 1 is  equivalent  to  v^  = 727  fps  (222  m/sec ) . The  predic- 
tions of  the  three  penetration  theories  are  compared  with  experimental 
3 9 

results  ® in  Figure  13.  For  the  experimental  results,  the  range  of 
solid  Reynolds  numbers  at  impact  is  1 5 R^  S 20  . The  difference  be- 
tween Goodier 's  results  and  those  of  the  shallow  and  deep  penetration 
theories  is  due  primarily  to  the  fact  that  Goodier  uses  Meyer’s  law  to 
characterize  target  resistance  during  the  projectile  embedding  process 
(References  3 and  5).  In  the  present  work,  Meyer’s  law  is  replaced  by 
an  equation  of  motion  which  is  based  upon  the  cavity  expansion  theory 
(see  Part  III,  paragraph  22),  resulting  in  reduced  penetration  depths  at 
the  lower  velocities. 


Shallow  Penetration  of  Concrete 


i^5.  The  penetration  of  concrete  has  been  the  object  of  extensive 
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Figure  13.  Nondimensional  final  penetration 
depth  versus  impact  velocity  for  steel  spheres 
and  202ii-T3  aluminum  target 


experimental  study  for  several  decades,  and  a number  of  empirical  pene- 
tration equations  have  been  developed  from  the  r^siilting  data.  For 
example,  a fairly  reliable  equation  for  the  penetration  of  reinforced 
concrete  by  various  types  of  projectiles  is  given  by^^ 


222  W D 


0.215 


^f  = 


A VT 


+ 2 ± 1% 


(Il8) 


where  ^ is  specified  in  inches,  W in  poimds,  D in  inches,  v^  in 

fps,  A in  square  inches,  and  Y in  psi.  No  range  of  applicability  is 

defined  for  this  equation,  but  sample  calculations  given  in  Reference  10 

lie  in  the  range  152  ra/sec  ^ v^  ^ 7^2  m/sec  (i.e.,  500  fps  ^ ^ 2500 

fps).  For  comparison  with  penetration  theory  predictions,  the  above 

equation  will  be  used  to  make  calculations  in  the  velocity  range 

100  m/ sec  < v < 1000  m/ sec  ( 328  fps  S v S 3280  fps ) , corresponding  to 
o ° 11 

0-7  5 S 70  at  impact.  The  concrete  static  properties  are 
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Pq  “ 150  pcf  (2.^0  ^/cn,'^) 

I - 5000  psi  (0.3^5  kliar) 
f. 

E — 4.10  X 10  psi  (283  tear) 


t 


12 


The  presswe- density  relation  illustrated  in  Figure  l4  represents 


Figure  l4.  Typical  pressure-density  relation 

for  concrete 


typical  'behavior  of  concrete  under  pressure.  The  static  value  of  Y 
is  used  in  tiie  theoretical  calculations  as  well  as  in  Equation  llS^  and 
the  resulting  theoretical  predictions  of  final  penetration  depth  then 
correspond  to  upper-hoxand  values.  However 5 if  a aiaximum  dynausie  value 
of  Y (i,e.5  40  percent  higher  than  the  static  value)  is  chosen ® the 
theoretical  penetration  depths  are  reduced  'by  only  about  15  percent. 
Since  Equation  II8  does  not  account  for  projectile  geometry^  two  sep- 
arate ogive  nose  shapes  are  used  in  the  theoretical  calculations , 

CRH  - 1.25  (L/D  - 1)  and  CRH  - 9.25  (L/D  = 3)  . 0th  erwises,  the  pro- 
jectile characteristics  are  chosen  to  be  W - 1000  lb  (4-55  kg)  and 
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D = 12  in.  (30.5  cm)  g -whicli  are  typical  of  the  projectiles  for  which 
Equation  l.l8  was  deduced.  Results  are  compared  in  Figure  15 » which 
shows  good  agreement  between  empirical  and  theoretical  predictions. 


J— 

u 


X 

s- 

0. 

UJ 

D 

z 

o 

sr 

h- 

UJ 

z 


UJ 

CL 

< 

z 


lOh- 


sh 


0^ 


0 ”*loo©  ' aooo  3000 

IMPACT  VELOCITY  j FT/SEC 


Figure  15.  Final  penetration  depth  versus  impact 
velocity  for  5000-*psi  concrete 


Perforation  of  an  Aluminum  Slab 


hG.  Wilkins^^  has  performed  2D  finite-difference  calculations 
analyzing  the  perforation  of  an  aluminum  slab  by  a sharp  steel  pro- 
jectile. This  sort  of  problem  is  of  particular  interest  in  the  present 
investigation  since  it  represents  an  extreme  example  of  penetration  in 
a layered  target.  The  first  layer  is  the  slab,  which  has  a finite  thick- 
ness as  well  as  a finite  density,  yield  strength,  and  elastic  modulus; 
the  second  layer  is  air,  which  has  a negligible  density  and  no  resis- 
tance at  all  to  shear  (i.e.,  the  second  layer  is  a fluid,  in  which 
= 00  for  all  q > O).  Wilkins  makes  two  separate  ealcula^ions ; 

(a)  in  his  first  calculation,  the  ultimate  strength  of  606I-T6  aluminum 
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Is  used  as  the  target  yield  strength;  and  (b)  in  his  second  calctilatioag 
this  yield  strength  is  doubled.  The  predictions  of  the  shallov  and  deep 
penetration  theories  for  layered  targets  will  now  be  compared  with 
Wilkins®  results.  TOie  aluminum  slab  has  a thickness  of  1 in.  (2.5^  em), 
and  the  rominal  target  properties  are 

g 

Pq  = 169  pcf  (2.T  gm/cm  ) 

Y = U3,500  psi  (3  kbar) 

E — 10*^  psi  (689  kbar) 

E*  - 0 
t 

and 

PqC^  — 10*^  psi  {689  kbar) 

Target  compressibility  is  adequately  characterized  by  Equation  for 
the  range  of  dynamic  pressxires  encountered.  The  projectile  characteris- 
tics are 


and 


W = 0.0183  lb  {8.32  gm) 
D = 0.3  in.  {O.T62  cm) 
L/D  - 1 (conical  nose) 

= 275^  fps  (8i^0  m/sec) 


U7.  Since  Wilkins  does  not  give  deceleration  histories  for  the 
projectile,  the  most  stringent  comparison  of  2D  and  theoretical  results 
lies  in  the  examination  of  residual  projectile  kinetic-  energy  (i.e., 
the  kinetic  ener^r  of  the  projectile  after  it  has  passed  completely 
through  the  slab).  Theoretical  and  2D  predictions  of  residual  kinetic 
energy  are  compared  in  nondimens ional  form  in  Figure  16,  where  is 

the  residual  projectile  velocity  and  v^^  chosen  as  the  non- 

dimensional  impact  velocity  (note  that  ~ at  impact ) . 

The  variation  of  nondimensional  impact  velocity  is  achieved  by  varying 
Y instead  of  v^  , and  the  predictions  of  kinetic  energy  loss  exceed 
Wilkins’  values  by  as  much  as  75  percent,  which  may  be  considered  "good” 


5^ 


Figure  l6.  Nondimensional  residual  kinetic 
energy  versus  impact  velocity  for  conical 
nosed  steel  projectile  after  perforating 
1-in. -thick  aluminum  slab 

or  "bad"  depending  on  one's  viewpoint.  (Wilkins’  results  for  Y = 3 kbar 

X3 

are  in  agreement  with  experiment.  ) This  example  indicates  to  some 

extent  the  accuracy  with  which  the  two  penetration  theories  predict 

projectile  deceleration  and  energy  loss  at  low  to  moderate  values  of 

near  an  interface  between  two  layers  which  have  extremely  different 

material  properties.  At  low  values  of  R , the  deep  penetration  theory 

s 

predicts  higher  energy  loss  than  the  shallow  penetration  theory,  but 

this  trend  is  reversed  as  R increases. 

s 

Deep  Penetration  of  a Multilayer  Earth  Target 

^8.  In  July  of  19T^,  the  Defense  Nuclear  Agency  (DNA)  sponsored 
a series  of  large-scale  projectile  penetration  tests  which  were  con- 
ducted by  Sandia  Laboratories  at  the  Watching  Hill  Blast  Range  near 

lU 

Ralston,  Alberta,  Canada.  The  three  test  projectiles  ranged  in  weight 
from  200  to  UOO  lb  (91  to  l82  kg),  final  penetration  depths  ranged  from 
30  to  99  ft  (9  to  30  m),  and  impact  velocities  ranged  from  306  to 
619  fps  (93  to  189  m/sec).  This  represents  one  of  the  best  docimiented 


55 


{ 


large-scale  earth  penetration  experiments  ever  conducted ® consisting 

of  eight  successful  shots  into  a target  which  had  “been  previously  sub- 
jected to  an  extensive  material  property  investigation.  A s\’mmaiy  of 
the  test  results  is  given  iu  Table  2,  which  is  taken  directly  from 
Reference  llj-. 

49.  The  shallow  and  deep  penetration  theories  will  now  be  used 
to  calculate  projectile  motion  in  the  idealized  four-layer  target  site. 
Bilinear  shear  properties  and  initial  densities  are  tabulated  for  the 
target  in  Table  3 (E  =0  for  all  layers).  These  properties  were  used 

^ 15 

in  a pretest  calculation  by  Rohani  and  are  obtained  from  the  results 

l4 

of  xinconfined  djTiamic  tests.  The  relation  between  density  and  pres- 
l4 

sure  is  given  for  each  layer  in  Figure  17 . The  projectile  weights 
and  diameters  are  200  lb  (9I  kg),  28O  lb  (12T  kg),  UOO  lb  (I82  kg),  and 
4*125  in.  (10.48  cm),  6 in.  (15.24  cm),  6.5  in.  (I6.51  cm),  respectively. 
The  400“lb  projectile  has  a 9-25  CRH  ogive  nose  shape  (L/D  = 3);  the 
other  two  have  6.0  CRH  ogive  nose  shapes  (L/D  = 2.4).  Comparisons  of 
theoretical  and  experimental  final  penetration  depths  are  given  in 
Figures  18-20,  and  a comparison  of  theoretical  and  experimental  decelera- 
tion records  for  Test  No.  6 is  shown  in  Figure  21.  These  calculations 

2 

fall  in  the  range  of  moderate  to  high  solid  Reynolds  number  (R  ~ 10  ) , 

s 

and  the  shallow  penetration  theory  exhibits  an  increasing  tendency  to 

underpredict  final  penetration  depth  at  the  higher  impact  velocities 

(i.e,,  at  the  higher  values  of  R ) . On  the  other  hand,  the  deep  pene- 

s 

tration  theory  reproduces  the  observed  final  depths  with  fair  accuracy. 
Examination  of  the  deceleration  records  in  Figure  21  reveals  that 
neither  theory  duplicates  the  details  of  the  experimental  deceleration 
curve  I however,  the  deep  penetration  theory  predicts  decelerations  which 
are  more  nearly  in  agreement  with  the  experimental  results. 


of  a Multilayer  Earth  Target 


50.  Some  of  the  deepest  earth  penetration  events  on  record  were 
observed  in  tests  conducted  by  Sandia  Laboratories  at  the  Tonopah  Test 
Range  (TTR)  near  Tonopah,  Nevada.  The  results  of  these  tests  are 
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Owantities  in  parentheses  were  derived  from  data  recorded  on  a high-frequency  early  time  recording  system.  The  other  decelera- 
tions tabulated  are  based  on  data  from  a low-frequency  on-board  digital  recorder. 


Table  3 

Material  Properties  for  Watchinf^  Hill  Test  Site 


Layer 

^0 

Y 

E 

Depth  of 
Layer  Bottom 

Mo. 

pef 

bar 

bar 

ft 

1 

93.0 

3.10 

276.0 

8 

88.6 

2.07 

138.0 

16 

<— i 
0 

116.0 

2.07 

138.0 

2h 

h 

123.0 

2.07 

82.7 

CO 

Figure  17.  Pressure-density  relations 
for  Watching  Hill  test  site 
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FINAL  PENETRATION  DEPTH,  FT 


DECELE 


I Vq»485  ft/sec  I 
•END  OF  RECORD 


15 

DEPTH, M 


DEPTH,  FT 


Figure  21.  Projectile  deceleration  versus  depth 
for  Watching  Hill  Test  No.  6 


reported  by  Patterson' 
characteristics ; 


for  a projectile  with  the  following 


W = 650  lb  (295  kg) 

D - 9 in.  (22.9  cm) 

CRH  - 6,0  (ogive  nose) 

IT 

Peterson  and  Hadala  have  made  upper-  and  lower-bound  strength  estimates 
for  the  Antelope  Dry  Lake  where  the  tests  were  conducted.  These  esti- 
mates were  made  using  available  unconfined  strength  data  and  visual 
classifications  for  the  Antelope  Lake  and  Main  Lake  beds  along  with 
constitutive  property  analyses  for  similar  playa  deposits  at  ITR  and 
at  the  Nevada  Test  Site.  The  estimated  strengths  are  listed  for  the 
idealized  four-layer  target  in  Table  U. 


Table  k 

Material  Properties  for  TTR  Antelope  Dry  Lake  Site 


Layer 

PQ 

Lower-Bound  Y 

Upper-Bound  Y 

Depth  of 
Layer  Bottom 

No, 

pcf 

bar 

bar 

ft 

1 

103 

0,689 

3,i*H8 

25 

2 

99 

0.001 

3M8 

37 

3 

111 

5,172 

10.31+5 

109 

h 

103 

3.Hlt8 

IO.3I+5 

GO 

For  the  upper-bound  strengths,  E 15O  Y ; and  for  the  lower-bound 
strengths,  E **  75  Y.  The  target  is  assumed  to  be  incompressible. 

51.  Experimental  and  theoretical  results  for  final  penetration 
depth  are  compared  in  Figures  22  and  23,  with  the  shaded  area  repre- 
senting the  range  of  predicted  depths  resulting  from  the  estimated 
bounds  on  the  target  properties.  The  difference  in  the  predictions  of 
the  two  penetration  theories  is  dramatic;  the  shallow  penetration 
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Figure  22.  Range  of  final  penetration  depth  versus 
impact  velocity  for  TTR  Antelope  Dry  Lake,  as  pre- 
dicted hy  the  shallow  penetration  theory 
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Figure  23.  Range  of  final  penetration  depth  versus 
impact  velocity  for  TTR  Antelope  Dry  Lake,  as  pre- 
dicted by  the  deep  penetration  theory 
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theoiy  exhibits  a growing  tendency  to  underpredict  final  depth  with  in- 
creasing irapact  velocity 5 while  the  deep  penetration  theory  predicts  an 
approximately  linear  relation  between  final  depth  and  impact  velocity 
which  is  in  rough  agreement  with  the  experimental  results.  Figure  23 
also  serves  to  illustrate  the  prediction  spread  which  is  to  be  expected 
when  target  property  estimates  are  made  v/ithout  the  benefit  of  a detailed 
investigation  of  the  target  site  itself.  In  these  calculations,  the 

solid  Reynolds  number  at  a velocity  of  2500  ft/sec  is  on  the  order  of 
3 

10  , No  experimental  deceleration  records  are  available. 


PAHT  ,V:  CONCLUSIOHS  AUD  RECOMffiMDATIOMS 


52,  As  a result  of  efforts  made  under  the  present  research  pro- 
grams the  projectile  penetration  theory  originally  formtilated  hy 

3 

Goodier  has  been  extended  so  that  the  theory  is  no¥’  applicable  for 
deep  penetration  as  well  as  for  shallow  penetration  in  both  homogeneous 
and  layered  cohesive  targets.  The  initial  formulation  of  the  theory 
was  based  on  an  ^ hoc  analogy  with  the  dynamic  expansion  of  a spherical 
cavity  in  an  infinite  medium.  The  final  step  in  the  theoretical  de- 
velopment was  directly  empirical  in  that  a single  function  i|;(R  ) had 

s 

to  be  chosen  which  relates  dynamic  penetration  resistance  to  the  solid 
Reynolds  number  as  well  as  to  the  projectile  geometry.  This  was  ac- 
complished by  means  of  a trial  and  error  comparison  of  experimental  and 
theoretical  results,  hut  it  is  emphasized  that  the  expression  for 

which  was  chosen  (Equation  8U)  is  the  same  for  all  targets  and 
projectiles.  The  deep  penetration  theory  now  represents  a self- 
contained  tool  for  making  penetration  predictions  and  requires  only  the 
projectile  characteristics  and  the  target  constitutive  properties  as 
input . 

53.  The  deep  penetration  theory  is  most  applicable  for  cohesive 
targets,  such  as  metal,  rock,  and  concrete.  On  the  other  hand,  it  is 
least  applicable  for  granular  targets,  such  as  sand,  which  exhibit  a 
strong  dependence  between  shear  strength  and  confining  pressure. 

The  comparisons  of  theoretical  and  experimental  results  for 
cohesive  targets  in  Part  IV  suggest  reliability  at  solid  Reynolds 

3 

numbers  as  high  as  R 10  , at  least  for  nose  shapes  with  L/D  >_  1 . 

For  L/D  < 1 5 the  range  of  reliability  probably  decreases  with  de- 
creasing L/D  . When  L/D  ^ 1 , the  theory  results  in  an  approximately 
linear  relation  between  final  penetration  depth  and  impact  velocity, 
which  is  consistent  with  experimental  observation.  The  theory  also 
predicts  an  approximately  linear  relation  between  final  depth  and  pro- 
jectile frontal  loading  W/A  , as  well  as  a nonlinear  relation  between 
final  depth  and  projectile  nose  length, 

55.  In  those  cases  where  it  is  applicable,  the  deep  penetration 


6k 


theory  caii  generally  be  expected  to  predict  peak  deceleration  within 
about  a factor  of  two  and  final  penetration  depth  within  about  20  per- 
cent when  the  target  density  and  unconfined  strength  are  known  with 
reasonable  a.ccuracy.  Other  quantities,  such  as  total  duration  of  the 
penetration  event,  will  probably  be  accurate  within  about  a factor  of 
two.  Error  in  prediction  will,  of  co\irse,  increase  as  uncertainty  in 
target  properties  increases. 

56.  Future  efforts  should  be  directed  toward  further  comparison 
of  theoretical  and  experimental  resxilts  and  toward  the  development  of 
a penetration  theory  for  materials  which  are  governed  by  the  Mohr- 
Coulomb  failure  condition  (shear  strength  proportional  to  confining 
pressure).  The  Mohr -Coulomb  condition  is  more  generally  applicable  to 
soils  than  is  the  Von  Mises  condition  which  is  used  in  the  present 
theory. 
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APPENDIX  A:  NOTATION 


a, a, a Radial  position,  velocity,  and  acceleration,  respectively,  of 

spherical  cavity  vallt 

Ej  Projectile  radius  at  layer  interface  location 

IfeximuBi  projectile  radius  - D/2  (Figtares  3 and  U) 

a Projectile  radius  at  target  surface 

^ 2 2 
A Maximum  projectile  cross-section  area  = ira^  = ttD  /h 

"b  Radial  position  of  spherical  plastic  front  (Figure  l) 

Initial  radial  position  of  particle  located  at  position  h 

®1’®2  Cavity  expansion  inertial  coefficients  (Equations  UO  and  Ul)tt 

C Dilatational  wave  speed 

CRH  Ogive  caliber  radius  (Equation  6?) 

D Maximum  projectile  diameter  = 2a^  (Figures  3 and  U) 

E Modulus  of  elasticity  in  compression 

Strain-hardening  mod-ulus  in  compression 

f ,f  Fionctions  of  integration  (Equations  9i  and  11) 

f^  Projectile  nose  shape  factor  (Equations  62,  63,  65,  72,  and  78) 
f^^j  Effective  nose  shape  factor  near  a layer  interface  (Equation  9^) 

F Net  axial  resisting  force  exerted  on  projectile  by  target 

^ 2 2 
g Gravitational  acceleration  = 32.2  ft/sec  (9.8  m/sec  ) 

h Radial  position  of  interface  between  two  concentric  layers 
(Figiu-e  1) 

h Initial  radial  position  of  interface  between  two  concentric 
layers 

H Distance  from  base  of  projectile  nose  to  layer  interface  loca- 
tion (Figure  5) 

L Projectile  nose  length  (Figures  3 and  I4) 

m Summation  index 

M Projectile  mass 


t A dot  above  any  quantity  denotes  differentiation  with  respect  to 

time. 

tt  With  the  exception  of  and  , the  subscripts  1 and  2 denote 

first  and  second  layer  quantities,  respectively. 
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Compressive  nonnal  stress  at  cavity  surface  (Equations  7 and  30) 
and  at  projectile  frontal  surface  (Equation  48) 

Dynamic  pressure  at  cavity  sxzrface  (Equation  38)  and  at  pro- 
jectile frontal  siirface  (Equation  49) 

Compressive  normal  stress  at  cavity  surface  due  to  material 
"behavior  in  shear  (Equation  39) 

Volume-averaged  dynamic  pressure  in  plastic  region  (Equations 
42  and  10?) 

Projectile  penetration  depth,  velocity,  and  acceleration, 
respectively 

Final  penetration  depth 

Eulerian  radial  coordinate 

Initial  radial  position  of  particle  located  at  position  r 
Sho'wn  in  Figure  4 

Solid  Reynolds  niunher  (Equation  68) 

Effective  distance  from  "base  of  projectile  nose  to  plastic 
front  (Figure  5) 

Time 

Outward  radial  particle  velocity 

Components  of  Vp  , respectively,  normal  and  tangent  to  projec- 
tile frontal  surface  (Figure  3) 

Projectile  impact  velocity 

Target  particle  velocity  adjacent  to  projectile  frontal  surface 
(Figure  3) 

Residual  projectile  velocity 
Projectile  weight 

Kondimensional  cone  position  (Equation  52) 

Equivalent  nondimens ional  cone  position  for  ogive  (Equation  55) 
Yield  strength  in  compression 

Axial  distance  measured  aft  of  projectile  nose  tip  (Figures  3 
and  4) 

Material  compaction  coefficient  (Equation  12) 

Kondimensional  material  parameters  (Equations  17  and  24) 

Discontinuity  in  a at  r = h(t) 

r 

Ogive  parameter  (Equation  67) 

Locking  strain  (Equation  l) 


A2 


*€  0) 


e Hcrmal  strains  in  radial  and  circumferential  directions  5 

respectively 

n Showi  in  Figxire  it 

C Hondimenslonal  depth  (Equation  115) 

TT  3.litl6 

p Material  density 

p n Locked  plastic  density 

p Projectile  density 

P 

Pq  Initial  density 

Limiting  value  of  as  r-^b  in  the  elastic  region 

^ Limiting  value  of  as  r->b  in  the  locked  plastic  region 

Normal  stresses  in  radial  and  circumferential  directions, 
respectively 

<|)  Cone  half-angle  (Figure  3) 

Equivalent  cone  half-angle  for  ogive 
Function  relating  v and  R (Equations  69  and  8U) 
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